Avoiding BBN Constraints on Mirror Models for Sterile Neutrinos 



R.N. Mohapatra and S. Nasri 

Department of Physics, University of Maryland, College Park, MD-20742, USA 

(Dated: July , 2004) 

We point out that in models that explain the LSND result for neutrino oscillation using the mirror 
neutrinos, the big bang nucleosynthesis constraint can be avoided by using the late time phase 
transition that only helps to mix the active and the sterile neutrinos. We discuss the astrophysical 
as well as cosmological implications of this proposal. 
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INTRODUCTION 

The existence of neutrino oscillations has now been 
confirmed for solar and atmospheric neutrinos as well as 
for reactor and accelerator neutrinos. It is remarkable 
that all the data from many different experiments can be 
well understood in terms of only three neutrinos that mix 
among themselves. They imply very narrow ranges of 
both the mass difference squares among these neutrinos 
as well as mixings. 

There is however another piece of evidence for os- 
cillations which if confirmed will require severe depar- 
ture from the successful three neutrino scheme just men- 
tioned. It is the apparent observation of the muon anti- 
neutrino oscillating to the electron type anti-neutrino in 
the Los Alamos LSND^J experiment. An attempt was 
made to confirm this result by KARMEN0 collaboration 
which eliminated a large fraction of the parameter space 
allowed by LSND. It is hoped that the Mini-BOONE ex- 
periment at FERMILAB currently under way will settle 
the issue in near future 0]. 

If the LSND experiment is confirmed, one straightfor- 
ward way to understand the results would be to postulate 
the existence of one or more extra neutrinos with mass in 
the eV range that do not interact with the W boson, the 
so-called sterile neutrinos and have them mix the known 
neutrinos. There have been various versions of this sug- 
gestion depending on the detailed mass arrangement of 
the sterile neutrinos with respect to the known ones: they 
are known in the literature as the 2+201, 3+1,0 as well 
as 3+20 models. Of the three, 3+1 model seems less 
disfavored than the 2+2 by the null results of other oscil- 
lation experiments. However the more recently proposed 
3+2 scenariojfj that involves two sterile neutrinos is ap- 
parently in better agreement wil all data than the others. 

The major challenge posed by the sterile neutrino for 
theory is to understand its ultra-lightness despite its be- 
ing a standard model singlet. A class of particle physics 
models that successfully answer this challenge are the 
mirror matter models. The basic assumption of these 
models is that there is an identical copy of the standard 
model (both constituents and forces) in nature0,0 that 
co-exists with the familiar standard model matter and 
forces. It is then clear that the same mechanism that 



keeps the active neutrinos light, will also keep the mirror 
neutrinos light which can therefore play the role of the 
sterile neutrinos. These models are inspired by the su- 
perstring theories and have been widely discussed 0, 0. 
Phcnomenological and astrophysical constraints on these 
models have also been extensively discussed [lol Hlf . 

Sterile neutrino models for LSND face two cosmologi- 
cal huddles that we would like to address in this paper. 
The issues are: how to make them consistent with (i) our 
understanding of big bang nucleosynthesis (BBN) and (ii) 
the recent bounds on neutrino masses from WMAP ob- 
servations. The first problem is that BBN allows the 
number neutrinos iV„, in equilibrium when the temper- 
ature of the Universe is one MeV is restricted by 4 He 
and £>2 observations to be very close to three|l2j. On 
the other hand for v s mass in the eV range and mixing 
in the few per cent range required to explain the LSND 
data, rapid v e — v s oscillations would lead to N v = 4 
for the 3+1 and 2+2 scenarios and N v = 5 for the 3+2 
scenario. 

The WMAP0 constraints are on the sum of all neu- 
trino masses in equilibrium at the epoch of structure for- 
mation which corresponds to a temperature around an 
eV. According to [lj[, J2 m v — 1-38 eV for one sterile 
and J2 m v — 2-12 eV for two extra ones assuming that 
they went into equilibrium at the BBN epoch. Thse con- 
straints are also quite important since taken at face value, 
they would seem to rule out the 3+2 model for LSND. 

It is therefore important to look for scenarios that may 
allow one to avoid both the above constraints while at 
the same time providing an explanation of the LSND ex- 
periment. Recently, it has been suggested that [la] by 
using late time phase transition to generate the masses 
and mixings of both the active and sterile neutrinos, one 
can avoid both these constraints. In ref.|l5l|. it is shown 
that this can be achieved by endowing two scalar fields 
(j> with vevs in the 100 keV range so that at the BBN 
time the sterile as well as the active neutrinos are mass- 
less. As a result there is no oscillation among them that 
can bring the sterile neutrinos into equilibrium. Since the 
sterile neutrinos decouple from Hubble expansion at very 
high temperatures, their abundance at the BBN epoch is 
suppressed leading to concordance with the BBN con- 
straints. Cosmological signatures of generic models of 
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this type have been given in Ref.[15j- 

In this paper we propose an alternative way to avoid 
the cosmological constraints using the same idea of late 
time phase transitions. We show that if the sterile neutri- 
nos are the mirror neutrinos, we need only generate the 
mixing between the active and the sterile neutrino (and 
not masses) by the late time phase transition to avoid 
the BBN and WMAP constraints. An advantage of this 
model is that the contribution of the sterile neutrinos to 
the energy density of the universe at the BBN epoch is 
giverned by a free parameter unlike the model of ref.|l5l|. 
We further find a convenient realization of mirror model 
with the seesaw scale in the TeV range which implies that 
we must employ the double seesaw mechanism |l6j to get 
small neutrino masses. We construct explicit scenarios 
with late phase transition and discuss their cosmological 
and astrophysical implications. The detailed field theo- 
retical models for them can be worked out but we do not 
discuss them here. 



AN EXTENDED MIRROR MODEL 

We start by reminding the reader about the generic 
features of the mirror models where one assumes that the 
universe consists not only of the observed standard model 
particles and forces but also coexisitng with an identical 
but different set of constituents experiencing analogous 
but different gauge forces. Gravity is however common 
to all the particles. The forces are dictated by the the 
gauge group G®G where one of the gauge groups G acts 
in the standard model sector and on its fermions whereas 
the other acts in the other and on the mirror fermions. 
The fermion spectra on both sides are identical. Mirror 
symmetry keeps the gauge couplings equal but the the 
effective strength of various forces in both sectors may be 
different due to different patterns of symmetry breaking. 
These models are inspired by the superstring theories as 
well as M-theory inspired brane-bulk models that have 
been widely discussed. 

As is clear, the neutrinos in the mirror sector do not ex- 
perience the known weak interactions and will not there- 
fore appear in the Z and W-decays. They can therefore 
play the role of the sterile neutrinos used in the inter- 
pretation of the LSND results. The mixing between the 
active and mirror neutrinos can arise in a manner con- 
sistent with gauge invariance (see below for details). For 
purposes of notation, we denote all particles and param- 
eters of the mirror sector will by a prime over the corre- 
sponding familiar sector symbol- e.g. mirror quarks are 
u', d', s', etc and mirror Higgs field as H' u d etc. 

Before proceeding further, let us discuss the origin 
of the masses and mixings for the active and sterile 
neutrinos. For this purpose, we extend the standard 
model gauge group in each sector to SU(2)l x U(1)i 3R x 
U(1)b-l, which is an anomaly free gauge group in 



the presence of the right handed neutrino v c . All the 
fermions have obvious quantum numbers under the gauge 
group. We add a gauge singlet chiral fermion, S in each 
sector, one per family. Mirror symmetry requires that 
we do the same in the mirror sector. In order to get the 
standard model gauge group from the extended group in 
each sector, we need to add a pair of new Higgs bosons 
A(l,+i,— 1) and a conjugate field A(l,— i,+l) in the 
visible sector and two similar fields in the mirror sector. 
We then add a gauge singlet Higgs field x ( and a mirror 
x'), which gives Majorana mass for the singlet fermions 
S by an interaction of the form \' a p{S a Spx + S' a S'gx'- 

The full superpotential relevant for neutrinos in each 
sector can be written as: 

W = KLH u v c + An/ 1 AS + A' SS X ; (1) 

where we have omitted an identical set of terms for the 
mirror sector and have suppressed the generation index. 
For three generation case that we will be interested in, 
h u , Ai and A' are 3x3 matrices. 

The U(l)i 3R x U(1)b-l part of the gauge symmetry 
is broken by the vev of field < A > assumed to be in 
the multi-TeV range. We choose the vev of x field to be 
in the GeV range. It is then clear that this leads to the 
double seesaw form 16] for the (v, i/ c , S) mass matrix: 

/ h v v \ 
M v =\ h T u v \ lVR . (2) 
V Af v R X' < x > / 

There is also a similar matrix for the mirror neutrinos. 
This leads to the light neutrino mass matrix of the form: 

M v = h v M- x X' < X > M R lT h T v v 2 (3) 

where Mr = Ai < A >. It follows that if we choose 
< x > about a GeV, h u — 1CT 1 and A' ~ 1CT 4 , then for 
Try- ~ tit- ~ 10~ 2 , the neutrino masses are in the 0.1 eV 

Mr M r 

range as required by observations. Also, typical neutrino 
mass textures can be built into the coupling matrix A'. 

A similar situation will occur in the mirror sector, 
where we can choose < %' > about a factor of 10 higher 
to get m Vs in the eV range to fit LSND (henceforth, we 
will call the sterile neutrinos v' as v s ). 

In order to generate mixing between the active and 
sterile neutrinos, we postulate the existence of a scalar 
field </> that mixes the two sectors. This can only be 
done through an interaction of the form f3SS'<j> [ljfl. A 
simple tree level diagram via the exchange of v c , S and 
v c ' ', S' then leads to an effective coupling of the form: 
A" L Af^A.r"^ ' wnere ^ — fth v ~ 10~ 2 . We assume as 
in ref.flijl that the vev of the field < <j) >~ 100 keV, so 
that at the BBN epoch the active and sterile neutrinos 
are unmixed 1 . The resulting v — v s mixing is then given 



The smallness of the <j> vev can be justified if we embed our theory 
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by: m v - Vs ~ A" ' ~ lO" 6 < > for A" ~ 10" 2 . 
This gives the right order of mixing for the LSND ex- 
periment. One difference between our model and that of 
ref.^| is that, late phase transition was used in 15] to 
generate both masses and mixings whereas in our case, 
only the mixing need to be generated at a late stage. 
As we will see, the mirror model has the advantage that 
contribution of the sterile neutrinos to the energy den- 
sity of the Universe at the BBN epoch is given by an 
arbitrary parameter. The model can therefore work even 
if the BBN constraints on the number of extra neutrino 
species tightened further. 



BBN, ASYMMETRIC INFLATION AND 
NEUTRINO MIXINGS 

Before we address the issue of neutrino mixings and 
BBN, we note that in the mirror model, we have three 
light neutrinos, a mirror photon and a mirror electron 
that could be potential contributors to the energy den- 
sity at BBN epoch and affect the success of BBN. In order 
to reduce their contribution to a negligible level, the idea 
of asymmetric inflation was proposed in the second paper 
of ref. , according to which it is assumed that the reheat 
temperature after inflation in the mirror sector is lower 
than that in the visible sector by a factor of 10 or so i.e. 
T' R ~ Tr/10. If the interactions linking the two sectors 
are such that they are not in thermal contact for T < Tr, 
then Hubble expansion will roughly maintain the ratio of 
the two temperatures till the BBN epoch apart from mi- 
nor corrections arising from particle annihilation in both 
sectors. Thus at the BBN epoch, the total contribution 
of the light mirror particles to ptot is at the level of about 
10 p 7 which therefore keeps the predictions of standard 
BBN unchanged. This will also keep the number densi- 
ties of the light particles such as mirror neutrinos and 
mirror photons suppressed. 



CONSTRAINTS AND CONSEQUENCES 

The first constraints on the parameters of the model 
come from the fact that at temperatures below the re- 
heating temperature after inflation, the interaction rates 
for L + H — > V + H' + <\> should be out of equilibrium, 
otherwise, the two sectors will be in the same thermal 
bath and at BBN, the effective N v will far exceed the 
allowed limit due to the fact that the population of the 
sterile neutrinos will build up their density to the level 



into a brane bulk scenario and have the < <j> > vev occur in a 
distant brane and get transmitted to the our brane via a bulk 
scalar field. 



of ordinary active neutrinos. Only exception is if the re- 
heat temperature after inflation is below an MeV which 
we do not invoke here. This translates into the following 
constraints in the three temperature regimes discussed 
below. 

(i) T >< H > , < H' > : In this region, the condition for 
being out-of-equilibrium is 

This inequality implies that for T 3> TeV (note that the 
mirror Higgs mass is expected to be in the TeV range), 
the visible and the mirror sector in our model are in equi- 
librium. So to be consistent with BBN requirements, we 
must require that the reheat temperature after inflation 
be less than about a TeV or so. 

(ii) 0.1 MeV <T << H >, < H' >: For T < M H < ~ 1 
TeV, the Higgs fields decay instantly In this regime, the 
effective interaction connecting the visible to the mir- 
ror sector is the coupling vv s 4> with a strength given by 

SW a ~ A" (-j^jp-j ~ 10~ 6 . In discussing whether 
this interaction is in equilibrium, it has been noted in 
ref.|l5l] that the process vv s — > 0, vanishes in the limit of 
rricj, = by energy momentum conservation. The rate 
for this process must therefore be proportional to m|(T). 
This effective thermal mass is given by n 2 T 2 /16tt 2 . If we 
choose the scalar self coupling n to be of order 10~ 2 , due 
to the fact that the number density of v s is down by a 
factor of 10 -3 , we expect the rate T(yv s — > <f) ~ 10~ 21 T. 
This leads to the conclusion that above T ~ 100 keV, the 
interactions connecting the visible with the mirror sector 
are out of equilibrium. Thus (i) and (ii) together then 
help to satisfy the BBN constraints. 

(iii) T < 0.1 MeV: 

This regime is below the scale of < <f> >. Therefore, 
v — v s — > (j) is kinematically forbidden since Imcj) becomes 
a pseudo-Goldstone boson. The only process that can 
lead to production of sterile neutrinos is vv — > v s v s . The 
rate for this process is given by 10~ 25 T. This process is in 
equilibrium below T < 100 eV. Below this temperature 
the mirror sector and the standard model neutrinos will 
thermalize without a significant transfer of energy. When 
the temeparture drops to T ~ m s , the sterile neutrinos 
decay into one of the three active neutrinos and the sin- 
glet Higgs and finally, when below the temperature ~ 1 
eV which is the mass of the <fi, it will annihilate via the 
process <fxf> — > vv leaving only a bath of active neutrino. 

To calculate the final temperature of neutrino bath in 
terms of the photon temperature, we first remember that 
the boson <j) is part of a supersymmetric multiplet whose 
scalar field part has only one surviving imaginary part 
and the real part has decoupled at very high temperature. 
We assume that the imaginary part is also superlight. 
Then we procedd through the following steps: 
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Just below T ~ m e , electron-positron annihila- 
tion heats up the photons leading to the relation 

T* = (^ij-) 1 ^ 3 ^. The u s and 4> at this stage are not 
in thermal contact with the active neutrinos. Once the 
temperature of the universe cools below 100 eV, u — u s <f> 
system comes into full thermal equilibrium. Using en- 
ergy conservation and taking the effect of the incomplete 
<j) supermultiplet into acccount, we get 



(3 
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(5) 



where n' is the number of sterile neutrinos and the factor 
take into account the contribution of the fermionic 
part of singlet Higgs superfield (f>. As the universe cools 
below the mass of u s , the u s decay to u +(f>. Using entropy 
conservation at this stage, we get 



11, 

y 



(3 + ^)1*., = (3 + ^+n')^ +0+l/s «n 



Using the above two equations, we get for the tempera- 
ture of the u + <j) system 



rpi _ 

u+ *~ 3+i 



Noting that p v oc (3 + t) T H>> 

we find the effective 
number of neutrinos at matter radiation equality is 



^-3[l + -^ Tr ] 1 /3 
T 



(8) 



When the temperature drops to T ~ the number 
of 0's get depleted via annihilation into standard model 
neutrinos. In this case the neutrino temperature is 



r 7 v ir L 



H/7,1/1 



/12 



(9) 



As an example if n' = 3 the effective number of neutrinos 
after BBN is N v = 4.08 and it is still constistent with 
CMB data. Future CMB experiments like PLANCK [T3| 
and CMBpol w ^ be able to improve the limit on N v 
and can provide a test of this model. The contribution 
of the active neutrinos to the critical energy density is 



53 v 1/3 m v 
2V 92h 2 



(10) 



Using the upper recent bound on the neutrino energy 
density on finds 



< 0.36 eV 



(11) 



We emphasize that this limit is only on the sum of the 
masses of the active neutrinos since in our model the 
sterile neutrinos have decayed away. 



Incidentally, the same steps can be repeated for theo- 
ries without supersymmetry. In which case, we will as- 
sume that the <p field has only a real part p and an imag- 
inary part x- The p field will have mass of order of the 
</>-vev or about 100 keV whereas we will assume the x 
mass to be an eV. In this case, below T ~ m e , first p 
decay dumps into the u + u s + x system giving 



1/3 



T, 



(12) 



1/4 



where T u+Us+<t> ~ ( 3+ | +n , ) T° due to u - u s - <p equi- 
librium. Taking entropy conservation at the subsequent 
u s decay and XX ~ * vv as m the previous case, we get 



1% 



1/3 



1/12 



(13) 



This changes the coefficient in the formula in Eq^J to 

1 /3 

(|y) and changes on the limit on the sum of active 
neutrino masses from 0.36 eV to 0.37 eV. 

Our scenario for sterile neutrinos has also interesting 
astrophysical implications. The first point is to look for 
any new mechanism for energy loss from the supernova 
core via emission of u s or <p. since u — u s mixing arises 
from spontaneous symmetry breaking at scale -C MeV, 
inside hot astrophysical environments such as a super- 
nova, the active and sterile neutrinos remain unmixed. 
As a result, there is no energy loss via the emission of u s . 
However, there could be energy loss due to the processes 
uu — > 4>(j),u s u s . The rates for these processes are esti- 
mated to be: ~ 10~ 25 T ~ (50 sec) -1 . Comparing this 
with typical supernova explosion time scale, we expect 
this energy loss mechanism not to be significant. Also 
due to zero mixing between u — u^ , all supernova results 
based on three active neutrinos |2lJ remain unaffected. 
Only in the very outer layers of the supernova explosion 
when the tempetature drops below 100 KeV, will these 
mixings become operative. 

DISCUSSION AND SUMMARY 

In summary, we have presented a mirror model for the 
sterile neutrinos that can explain the LSND results and 
yet be consistent with stringent constraints from big bang 
nucleosynthesis as well as cosmic microwave background 
as well as structure formation bounds on neutrino proper- 
ties. We make predictions for the effective neutrino num- 
ber to which the next generation CMB measurements are 
sensitive. An important requirement of this model is that 
the reheat temperature after inflation must be less than a 
TeV. The model has also other interesting properties dis- 
cussed earlier such as the mirror hydrogen being a dark 
matter 22] , which remain unaffected by our modification. 
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Similarly, suggestions that the ultra high energy neutri- 
nos could be originating from topological defects in the 
mirror sector |23j remain unchanged by our extension. 
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